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ABSTRACT 

We investigate the history of 6 Li and Li in population II dwarfs during the pre main sequence 
and main sequence. The evolution is followed using the CESAM code and taking into account the 
most recent physics. The effective temperature ranges from w 4700 K to « 6400 K and therefore 
concerns objects on the so-called Spite plateau and cooler. We find the 7 Li pre main sequence 
depletion is unable to account for the observations in the halo whatever the effective temperature. 
This supports microscopic diffusion and an additional non standard mixing process both acting 
during the main sequence. On the contrary the models 6 Li pre main sequence depletion appears 
too strong and is marginally compatible with recent detections. 

During the main sequence we introduce the effects of tachocline diffusion. This process is a 
rotationnally induced mixing acting at the top of the radiative core. We show that the differences 
in the early rotation history cannot result in scattered lithium abundances on the Spite plateau. 
Moreover the tachocline mixing process predicts 7 Li abundances in good agreement with the 
observations. We briefly address the question of turn off 7 Li poor stars. A modest accretion of 
lithium free matter would be enough to explain their low abundance in this element. 

We expect the scatter in 7 Li abundances is correlated to variations in [Fe/O] ratio for dwarfs 
cooler than 5500 K. Finally the tachocline mixing is robust with respect to the recent 6 Li obser- 
vations around the turn-off. We similarly suggest the [Fe/O] should be higher in objects with 6 Li 
and effective temperature below 6000 K. 

Subject headings: stars: population II — stars: light elements — stars : internal mixing 

on their possible surface abundance changes since 
they were formed. In particular 7 Li has been ex- 
tensively studied in these objects both from theo- 
retical and observational viewpoints. 7 Li is one of 
the rare nucleids that was synthetised through Big 
Bang Nucleosynthesis (hereafter BBN) in an ob- 
servable amount. Estimating its initial abundance 
from observations of the oldest stars provides es- 
sential indications on the conditions in the primor- 
dial Universe. More than 20 years ago Spite & 
Spite (1982) found that most halo field stars with 
[Fe/H] < -1.5dex and 5500K < T cff < 6300K dis- 
play similar 7 Li abundances. Since then this result 
was confirmed by many other observations (Hobbs 
& Duncan 1987; Spite & Spite 1993; Thorburn 
1994; Ryan et al. 1996). The 7 Li abundance on 
the so-called Spite plateau is now estimated to lie 
around 2.0-2.2 dex with a very small intrinsic scat- 
ter below 0.1 dex. The key issue for stellar physi- 



1. Introduction 

Lithium, beryllium, and boron isotopes get de- 
stroyed by proton capture in stellar interiors be- 
tween 2 10 6 K and 5 10 6 K. Therefore their sur- 
face abundances allow direct measurement of the 
depth of the outer convection zones. In addition 
when the temperature near the base of a convec- 
tion zone (hereafter BCZ) is too low they permit 
to set constraints on any process that would mix 
the matter beneath. The light elements have been 
extensively observed in main sequence low mass 
stars. With asteroseismology they are currently 
the only means to constrain low mass dwarfs' in- 
teriors. 

Population II main sequence stars being the 
oldest non-evolved objects, they should reflect the 
initial chemical conditions of our Galaxy. It is 
therefore particularly interesting to set constraints 



cists is the determination of how this abundance 
has changed since the halo stars were formed. 

Recently, new spectrographs such as UVES on 
the VLT have considerably increased the number 
of measurements in other light elements isotopes 
such as 6 Li (Asplund ct al 2005a, 2005b) and 9 Bc. 
This makes the determination of their evolution 
from the modelling point of view relevant. In con- 
trast to 7 Li, these nucleids are not produced in a 
significant manner through BBN but are the result 
of fusion or spallation processes in the interstel- 
lar medium as suggested by the pioneering work 
of Reeves et al. (1970). Beside this they have 
different depletion temperatures than 7 Li so that 
they provide complementary information on stel- 
lar structure and evolution. 6 Li is especially in- 
triguing as present measurements suggest it does 
not evolve over a broad range of metallicity (As- 
plund et al. 2005a) whereas it is expected to in- 
crease substantialy due to production by low en- 
ergy cosmic rays (Vangioni-Flam et al. 1999). 

In this work we present the predictions of mod- 
els of population II low mass stars on the 6 Li and 
7 Li surface abundances. The models are evolved 
up to 13 Gyr, a plausible age for halo stars hav- 
ing [Fe/H]=-2 dcx (VandenBerg 2000). The his- 
tory of the isotopes is investigated both on the pre 
main sequence (hereafter pre MS) and the main 
sequence (hereafter MS) using the most recent 
physics in terms of nuclear reaction rates, equation 
of state and opacities. We have, moreover, studied 
the effects of the rotationally induced tachoclinc 
mixing. This process has given encouraging results 
for the evolution of light elements on the surface 
of population I solar analogs but has never been 
studied in population II objects. The paper is or- 
ganized as follows : in §2 we describe in detail the 
general inputs to our code and the assumed ini- 
tial composition. The rotation history has been 
followed in detail. In §3 we present our assump- 
tions on the rotation history and consequently on 
the tachocline mixing process. Sections 4 and 5 
respectively investigate the effects of pre MS and 
MS history on 6 Li and 7 Li. Finally we discuss our 
results in §6. 



2. The evolutionary code and general in- 
puts 

We use the CESAM code to perform the com- 
putations (Morel 1997). This hydrostatic oncd- 
imensional code has been extensively employed 
to model various stellar types and stages of evo- 
lution over the last decade. It has also been 
used in solar structure computations and helio- 
/asteroseismology 1 . Various adaptations have 
been made to the code in order to use the most 
adequate inputs for the present issue. 

The equation of state: in the case of population 
II objects the metals' impact on the equation of 
state (EOS) is negligible. Consequently we use the 
OPAL2001 EOS tables for pure hydrogen-helium 
mixtures 2 (Rogers & Nayfonov 2002). 

The opacities: consistent with the EOS, we use 
the OPAL opacities. The opacity table consid- 
ered is the a-element enhanced table of F. Allard. 
For this table [a/Fe] = 0.3dex where a stands for 
the oj-elements (see hereafter). Low temperature 
opacity tables are based on computations similar 
to Alexander and Ferguson (1994). Below log T 
=3.75 the OPAL opacity tables are replaced by the 
tables provided by J. Ferguson (private communi- 
cation) . These tables were generated for the same 
intermetallic ratios as those of the OPAL opacity 
tables. 

The nuclear reactions: we use the NACRE 
compilation of nuclear reaction rates (Angulo 
et al. 1999). Our network includes the proton- 
proton chains and CNO cycle. Furthermore the 
6 Li(p, a) 3 He and 9 Bc(p, 2 H)2 4 He reactions are 
taken into account. 

The atmosphere modelling: the outer boundary 
conditions to the stellar structure equations are 
provided using the Nextgen atmosphere models. 
The atmospheric temperature-optical depth rela- 
tions have kindly been computed by P. Hauschildt 
for the required composition and on a grid that 
spans the effective temperature-surface gravity do- 
main we encounter. 

The convection: the convection zones are as- 
sumed to be fully homogeneous, and modeled us- 
ing the mixing-length theory (hereafter MLT) in 

lr The list of the scientific publications using CESAM is avail- 
able at http://www.obs-nice.fr/cesam/ 
2 Available at http://www-phys.llnl.gov/Rcscarch/OPAL/ 
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a formalism almost identical to that of Bohm- 
Vitense (1958). This formalism is precisely de- 
scribed in the appendix of Piau et al. (2005). We 
have considered «mlt = 1-766 which is our Sun- 
calibrated value. 

The diffusion processes: The microscopic diffu- 
sion is taken into account following the Michaud 
& Proffitt (1993) prescriptions. The details of the 
tachocline diffusion process are given in §3. 

In addition to these different points the ini- 
tial composition is an essential question. For in- 
stance, the impact of mctallicity on the extension 
of convection during the pre MS is a priori un- 
known. This impact is strong in the case of pop- 
ulation I object (Piau & Turck-Chieze 2002, here- 
after PTC02). Moreover in MS models changing 
the intermetallic ratios affects the temperature at 
the base of the convection zone for a given effec- 
tive temperature. Both phenomona could signif- 
icantly change the light elements' histories. We 
assume [Fe/H]=-2.0 dex and [a/Fc]=+0.3 dcx. a 
stands for the oxygen and the usual a-elements 
(Nc,Mg,Si,S,Ar,Ca and Ti). We consider a similar 
ratio for nickel and chromium as for iron while 
manganese and aluminum are supposed under- 
abundant by 0.15 dex and 0.30 dex with respect to 
iron (Gratton 1989; Magain 1989). Recent obser- 
vations of the [C/O] vs [O/H] function in halo stars 
(Akerman et al. 2004) suggest that the carbon is 
underabundant with respect to oxygen there. The 
deficiency being around 0.4 dex for the oxygen 
fraction accounted in our models it would bring 
the carbon to iron ratio close to solar. We have 
thus considered [C/Fe]=0 which is also in agree- 
ment with other studies regarding carbon in pop- 
ulation II (McWilliam 1995). 3 Nitrogen is as- 
sumed to follow the same trend as carbon so that 
[N/Fc]=0. The total helium mass fraction along 
with the ratios of 2 H, 3 He, and 7 Li to 1 H are set to 
their primordial values following the recent BBN 
calculations of Coc et al. (2004). We consider that 
6 Li=1.13 dex and 9 Be=-0.15 dex initially which is 
adapted from the observations of turn-off popula- 
tion II stars by Cayrel et al. (1999) and Pasquini 
et al. (2004) respectively. In both last cases we 

3 The situation seems different in extremly metal poor stars 
(-4.1 < [Fc/H] < -2.7dcx) where [C/Fc] or [(C+N)/Fe] 
is close to 0.2 dex (Cayrel et al. 2004, Spite et al. 2005). 
However, only the objects with [Fe/H] around -2.0 dex are 
within the scope of this work. 



increase the observed abundances by 0.2 dex in or- 
der to compensate for the microscopic diffusion at 
the halo age. 9 Be surface abundance is followed, 
however because of the higher temperature for 9 Be 
proton capture it is never affected by nuclear re- 
actions on our range of T e ff. We therefore shall 
not discuss 9 Be any further in this study. The 
tables 1 and 2 provide the detailed composition. 
Hydrogen and helium mass fractions are X= 0.75 18 
and Y=0.2479 respectively. Unless explicitly men- 
tioned we always adopt this composition. 

3. The rotation 

Population II stars cannot be observed during 
their evolution. Thus it is impossible to observe 
their surface rotation history. Currently all the 
main sequence population II stars exhibit very low 
equatorial velocities. The line broadening due to 
these rotation rates is comparable to the spectro- 
graphs resolution today and for most objects only 
upper limits of the projected equatorial velocities 
are accessible. Studying 9 halo stars Smith et al. 
(1998) found vsini < 3km.s~ 1 for all of them. 
Even this upper limit was lowered to 2 km/s by 
more recent observations (see Fig. 1 of Ryan et al. 
2002). Interestingly, the latter authors also report 
that 3 stars out of their sample exhibit slightly 
higher vsini in the 5.5 to 8.3 km.s -1 range. These 
peculiar stars are spectroscopic binaries. They ex- 
hibit 7 Li surface abundances at least 0.5 dex below 
the other objects and should not be considered 
when studying the isolated stars general lithium 
history we address here. However we briefly touch 
on the question of such Li-poor objects in §5. 

We compute the rotation evolution by applying 
different assumptions to the angular momentum 
losses and internal transport for a 0.7M Q model. 
With T cff = 5890K at 13 Gyr this star lies in the 
middle of the Li-plateau. Four different models 
of rotation were built. In the following subsec- 
tions we present our assumptions on the rotation 
history. First we discuss our choices of surface 
angular momentum loss laws that in turn deter- 
mine the surface rotation. Second we address the 
question of the inner radiation core rotation that 
is the relevant quantity for the tachocline diffu- 
sion efficiency Third we briefly describe the way 
tachocline mixing was implemented. 
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Table 1: Metal fractions relative to the Sun. 



Metals 


[X/H] 


0,Ne,Na,Mg,Si,P,S,Cl,Ar,Ca,Ti 


-1.7 


Fe,Ni,Cr,C,N,K 


-2.0 


Mn 


-2.15 


Al 


-2.30 



Table 2: Light element number abundances relative to hydrogen adapted from Coc et. al (2004) except for 
6 Li and 9 Be (see text). Following Coc et. al (2004) we consider the helium mass fraction to be Y=0.2479. 



ratios 






SLi/iH 


'Li/ 1 !! 


SBc/iH 


values 


2.60 10" 5 


1.04 10" 5 


1.34 10~ n 


4.15 10~ 10 


7.07 10~ 13 


values (log) 


-4.58 


-4.98 


1.13 


2.62 


-0.15 



3.1. The surface rotation 

We make the assumption that population II low 
mass stars experience an evolution of surface ro- 
tation which is similar to their population I coun- 
terparts: 

i) The star initially rotates as a solid body and 
following Bouvier et al. (1997) (hereafter BFA97) 
we consider an initial rotational period of P =8 
days. This period remains unchanged until the 
star decouples from its initial accretion disk which 
we assume occurs at Td =3 Myr (case A) or 0.5 
Myr (case B). 

ii) Once the star is no longer locked to its disk 
the angular momentum decreases following the 
Kawaler (1988) prescription. 

^ = -™ 3 <i|>' ,2 <^>- 1/2ifSi < si « (1 » 



5f = - K aa,J(-^yi*({L)-VHta > «,„ 

at Rq Mq 

(2) 

K is set to K Q = 3.25 10 47 g.s.cm 2 (cases A and 
B) which leads to the actual solar rotation for a 
solar model. For case C we take K = K /1O. 
This second value yields a 1.6km.s _1 equatorial 
velocity at 13 Gyr which is nearly the maximum 
velocity allowed by current observations (vsini < 



2km. s _1 ). The surface magnetic activity saturates 
around lOf^©. The actual threshold f2 sat = 14O 
(3.78 10 _5 rad.s _1 ) is accurately tuned from the ro- 
tation observations in middle aged open-clusters 
(BFA97). 

iii) To compute the rotational rate inside the 
radiation zone we explore two possibilities regard- 
ing its coupling to the convcctivc envelope. Firstly 
we consider that these regions are instantaneously 
coupled so that the whole star rotates as a solid 
always (cases A, B, C). Secondly we consider the 
radiative and convective zones to rotate as solids 
but not necessarily with the same speeds (case D, 
see the next subsection). 

The four empirical parameters (P ,Tdi O sat ,K) 
in the rotation law we adopt are directly calibrated 
from the observations on our Sun and popula- 
tion I solar analogs. It seems reasonable to as- 
sume similar angular momentum losses through 
magnetic braking for solar analogs and Li-plateau 
stars. The reason is that, in terms of surface condi- 
tions or convection zone extension, the Sun is quite 
comparable to our 0.7M Q population II [Fe/H]=- 
2 dex model. At 13 Gyr we find for this object 
: T off = 5890 K, surface gravity = 1.36 g Q and 
0.764 Ft,* as radius for the radiative core vs 0.713 
R for the actual Sun (global properties of plateau 
stars as a function of their masses can be found in 
Proffitt & Michaud 1991). We consequently ap- 
ply the braking law parameters from BFA97. The 
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Fig. 1. — Radiation zone angular velocity vs age 
for a O.7M model. The solid body rotation mod- 
els are shown by the solid line : K = K Q , and 
Td = 3 Myr (case A); the dotted line : K = K Q , 
and Td = 0.5 Myr (case B); the dashed line : 
K = K Q /10, and r d = 3 Myr (case C). Dot-dashed 
line represents the differential rotation model with 
K = K Q , r d = 3 Myr and r c = 1 Gyr (case D). The 
difference between the case A and case B models 
is related to early conditions and vanishes after ~ 
500 Myr. 

shorter P Q and/or T d the higher zero age main se- 
quence (hereafter ZAMS) rotation velocity (figure 
1). Then a higher f2 sat implies a slower pace for 
rotational slow-down during early MS. To study 
the effects of fast ZAMS rotation we built a model 
with Td = 0.5 Myr. We stress however that none 
of the three parameters (P , Td and O sat ) has a sig- 
nificant impact on rotation after ~ 0.5 Gyr where, 
in contrast, the K parameter plays a major role. 
K determines the star's asymptotic velocity. Con- 
sidering Td = 3 Myr, we explore its impact when 
it is varied from K Q to K Q /10 for which it causes 
surface rotation velocities to vary from ~ 0.6 to 
~ 1.7km.s _1 at 13 Gyr (figure 1). 

3.2. The internal rotation 

The rotational rate in the upper radiation zone 
is the relevant quantity for the tachoclinc diffu- 
sion computation. This rate is also the surface 
rotation rate when the rotation is solid. In the 
case of population I solar like stars there is an ac- 
cumulation of observational and theoretical clues 
suggesting that this should be correct during most 



of the main sequence. Helioseismic observations 
show that our Sun rotates nearly as a solid at 
least down to 0.2R Q (Thompson et al. 2003). 
As shown in Piau ct al. (2003), the surface ro- 
tation history in the young clusters does not seem 
compatible with differential rotation with the ra- 
dius. In this last work we came to the conclusion 
that solid body rotation was certainly achieved as 
early as the Hyades age. Furthermore one would 
expect more significant surface magnetic activity 
in young slow rotators if they had fast spinning 
cores (BFA97). Finally some recent models of an- 
gular momentum transport by internal waves sug- 
gest an efficient rigidification of the stellar interior 
both for the Sun (Talon, Kumar & Zahn 2002) and 
for the lithium plateau stars (Talon & Charbonnel 
2004). 

In the case of lithium plateau stars however 
we cannot completely exclude differential rotation 
with depth because part of the low mctallicity 
horizontal branch (hereafter HB) stars show sig- 
nificant surface rotation while on the other hand 
the turn off halo stars are extremely slow rota- 
tors. We therefore built a rotation model assum- 
ing K = K Q , Td = 3 Myr and t c = 1 Gyr. For 
this model, during the time interval dt, the radia- 
tive core and the convective envelope exchange an 
angular momentum fraction dJ = ^dt. 



Icorc H - Icnv 

is the angular momentum exchange that would 
achieve synchronisation of the core and the en- 
velope, t c is therefore a coupling timescale for ro- 
tation between these regions. As shown on fig- 
ure 1 this option yields to a significant increase 
of the rotation rate in the radiation zone. At 
13 Gyr the core of the 0.7 M Q star rotates at 
1.35 10 _5 rad.s _1 while the surface rotation veloc- 
ity linear is 1.5km.s _1 . 

For this model the radiative core rotates ~ 5 
times faster than the envelope. Indeed this high 
ratio could be expected from constraints on fast 
rotating HB stars. The coolest HB stars (T e ff < 
12000K) are spinning as fast as 40km. s _1 (Bchr 
2003). Sills and Pinsonneault (2000) made a sys- 
tematic study of the angular momentum history 
during the red giant branch evolution. Their mod- 
els save the greatest amount of angular momen- 
tum when specific angular momentum is constant 
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in the convective envelope and fully retained in 
the radiative core. For this situation the au- 
thors demonstrated that the fast rotating HB stars 
could be accounted for if turn-off progenitors ro- 
tate at 4km. s _1 . The subsequent observations of 
Ryan et al. (2002) suggest vsini < 2km/. _1 how- 
ever. Thus if these turn-off stars rotate as solid 
bodies their low angular momentum seems diffi- 
cult to reconcile with the high rotation velocities 
of the HB stars even for very favorable conditions 
of angular momentum distribution. We can es- 
timate the rotational velocity of the core for a 
0.8M Q turnoff star under the following hypothe- 
ses: i) The total angular momentum of the star is 
similar to its HB value : Jhb = Jcorc + Jonv Jhb, 
J core and J cnv respectively stand for the angular 
momenta of all the HB star, its turn off progeni- 
tor core and envelope, ii) Both the whole HB star 
and the radiative core of its turn off progenitor 
rotate as solids, iii) The surface velocities of the 
HB star and its turn off progenitor are respectively 
40km. s -1 and 2km.s~ 1 . 
Then 

__ IHB^HB — Icnv^cnv 
^corc j \3) 

Icore 

At 11.8 Gyr when central hydrogen exhaustion 
occurs (X c = 10" 3 ) I corc = 4.0810 53 g.cm 2 , I cnv = 
3.7110 51 g.cm 2 , and ft cnv = 2.45 lO^rad.s- 1 . 
The rotation rate of the core is tllUS ^corc — 
3.3810 _5 rad.s _1 . Now if we compute the rota- 
tional evolution of this 0.8M Q model with the 
case D parameters we find fi coro = 4.710~ 5 rad.s~ 1 
at the same age. Both evaluations are quite com- 
parable. 

Before moving to the rotation effects in the next 
sections let us sum up our four different assump- 
tions on rotation and what motivates them : 

Case A): Solid body rotation, K = K©, and 
Td = 3 Myr. This rotation model is calibrated 
thanks to the Sun and population I observations. 
Unless explicitly mentioned we will always use this 
rotational law. 

Case B): Solid body rotation, K = K Q , and 
r d = 0.5 Myr. 

Case C): Solid body rotation, K = K /1O, and 
Td = 3 Myr. Cases B) and C) evaluate the impact 
of increased rotation rates with respect to case A) 
either because of lower angular momentum losses 
or higher ZAMS rotation speed. 



Case D): Differential rotation model K = K©, 
and Td = 3 Myr. The coupling time scale between 
the convective envelope and the radiative core is 
t c = 1 Gyr. This model aims at evaluating the 
effects of a possible fast rotating core. 

3.3. The tachocline mixing 

The transition region between the solar convec- 
tive and radiative zones shows sharp variations of 
the rotation velocity with depth and has there- 
fore received the name of tachocline. There ro- 
tation and differential rotation with latitude in- 
duce a slow mixing (Spiegel & Zahn 1992). The 
tachocline mixing has been introduced in mod- 
elling studies about the Sun and population I 
solar-like stars (Brun ct al. 1999, 2002; Piau et 
al. 2003). It produces better agreement between 
theoretical and observed sound speed and explains 
the 7 Li and 9 Be main sequence histories. How- 
ever, the tachocline mixing effect has never been 
systematically studied in population II stars. 

We presented in detail our assumptions on the 
global rotation history hereabovc. The differential 
rotation f2 with latitude in the convection zone 
is then deduced assuming it is the solar differen- 
tial rotation for the solar rotation rate and oth- 
erwise scales as fi 07 (Donahue, Saar & Baliunas 
1996). It is noteworthy that the tachocline mixing 
has features making it robust with respect to the 
present constraints on stellar rotation: it is not 
related to differential rotation with depth which 
presently seems excluded in MS stars by observa- 
tions as well as by new models of internal angular 
momentum transport. Besides this and contrary 
to models where angular momentum is extracted 
by a diffusive process, it infers no 9 Be depletion 
in solar analogs. Our tachocline mixing prescrip- 
tion has two free parameters that must be empir- 
ically calibrated: the width of the tachocline and 
the buoyancy frequency in the tachocline region. 
The width is calibrated from solar seismic mea- 
surements to 2.5% of the total radius. The value 
of the buoyancy frequency is estimated to be 10 
/iHz. A change to 2 ^iHz is briefly investigated 
however (figure 4). A discussion of the choice of 
these values can be found in Piau et al. (2003). 

We do not take into account the tachocline mix- 
ing during the pre main sequence (which is con- 
sidered to end at 200 Myr) unless explicitly men- 
tioned. The reason is that the time for effective 
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diffusion induced by the tachociine process is the 
order of : 

tdiff = 450(-)" 1 Myr 

(Zahn 2004). Even in the case of extreme differ- 
ential rotation (fl = ft) this duration thus already 
exceeds 200 Myr. A detailed description of the 
tachociine process equations can be found in Zahn 
(2004). 

4. Lithium isotopes during the pre main 
sequence 

During the pre MS, low mass stars presumably 
evolve from a fully convective state to a radiative 
core and convective envelope structure. The tem- 
perature at the BCZ reaches a maximum which 
may induce light element depletion visible at the 
surface (see PTC02 and references therein). The 
lower the mass (i.e. the lower the ZAMS T e g), 
the deeper the convection zone at any stage of the 
evolution and the stronger the light element de- 
pletion. We evaluate the pre MS depletion of 6 Li 
and 7 Li based on models including no tachociine 
diffusion. The impact of the pre MS evolution 
is considered at 200 Myr. At this age the proton- 
proton chains provide more than 99% of the energy 
while the microscopic diffusion is negligible for all 
our models. We consider therefore 200 Myr as the 
end of pre MS. The computations were made for 
the helium mass fraction Y=0.2479 and the metal- 
licities [Fe/H]=-3 dex, [Fe/H]=-2 dex (standard 
case), [Fe/H]=-1 dex and [Fe/H]=-0.7 dex. As 
shown in table 3 and figure 2 no significant 7 Li de- 
pletion is induced during pre MS even for objects 
of ZAMS effective temperature at the cool end of 
the Spite plateau (~ 0.65M©). Below this limit we 
estimate a depletion phenomenon rapidly increas- 
ing as the effective temperature decreases. These 
predictions are however unable to fully explains 
the T c ff — 7 Li pattern presently observed below 
5500 K. Firstly in the T e ff — 7 Li plan, we compute 
a 0.18 dex 7 Li decrease per 100 K between 5175 
K (0.62M© model) and 4975 K (0.58M© model). 
This slope is smaller than the estimate of 0.27 dex 
based on observations between 5000 K and 5500 
K suggested by Ryan & Deliyannis (1998). Sec- 
ond the predicted depletion rate around 5000 K 
is clearly lower than what is observed. Finally 
even the current intensity computed for pre MS 
7 Li depletion could be questioned: in the case of 
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Fig. 2.— 7 Li pre MS depletion at 200 Myr in 
[Fe/H]=-2 dex (solid line) , [Fe/H]=-1 dex (dashed 
line), [Fe/H]=-3 dex (dash-dotted line) models. 
These models all exhibit similar pre MS 7 Li deple- 
tion patterns. The [Fe/H]=-0.7 dex models (dash 
three dotted line) exhibit higher depletion because 
of the increasing impact of metals on opacity. Ob- 
servations (diamonds) and upper limits (triangles) 
are from Ryan et al. (1996), Ryan & Deliyannis 
(1998), Ryan et al. (1999) and Ryan et al. (2001). 
For the sake of clarity error bars are not shown. 

the population I G-type stars it clearly exceeds 
what is observed (Ventura et al. 1998; PTC02). 
We conclude that the T ff — 7 Li relation currently 
observed in the halo stars and for temperatures 
cooler than the cool end of the Spite plateau was 
(at least partially) built after the pre MS. 

It is interesting to note that the computed pre 
MS 7 Li depletion does not depend on the metal- 
licity from [Fe/H]=-3 to [Fc/H]=-1 dex. As shown 
on figure 2 the predicted 7 Li — T e ff relations are al- 
most identical on this range of metal fraction. This 
point is a marked difference from the pre MS pop- 
ulation I depletion predictions. Contrary to pop- 
ulation I objects the metals represent a very small 
fraction of the plasma components in the halo 
dwarfs. Despite the numerous bound-bound and 
bound-free interactions with the radiation field the 
contribution of heavy chemical species to the opac- 
ity budget near the BCZ becomes small below 
[Fe/H]=-ldex. We computed the metal contribu- 
tion to the opacity in the [Fe/H]=-0.7, -1 and -2 
dex cases. For this purpose we used the monochro- 
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Table 3: ZAMS (age = 200 Myr)lithium to initial lithium ratios as function of temperature for microscopic 
diffusion models : [Fe/H]=-2 dex and Y=0.2479. 



Mass (M Q ) 


0.80 


0.75 


0.70 


0.65 


0.62 


0.58 


ZAMS T cff (K) 


6144 


5893 


5621 


5341 


5177 


4975 


7 LizAMs/ 7 Lio 


0.97 


0.96 


0.93 


0.83 


0.69 


0.30 


6 LizAMs/ 6 Lio 


0.43 


0.15 


1.1 IO" 2 


3.6 10" 6 


< io- 7 


< io- 7 



matic calculations of Iglesias & Rogers (1996). 
The contribution of a component is computed by 
subtracting the opacity evaluated for the mixture 
without the component to the opacity of the global 
mixture. This requires that the component whose 
opacity contribution is evaluated does not provide 
a significantly fraction of the free electrons 4 . Ta- 
ble 4 compares these computations showing the in- 
creasing role of metals in the opacity budget with 
increasing [Fe/H]. 

Being more fragile, 6 Li is more strongly affected 
by the pre MS peak temperatures near the BCZ. 
We compare here our computations to data col- 
lected in Nissen et al. (1999) and mostly Asplund 
et al. (2005a, 2005b). For this sample we restrict 
ourselves to the nine stars showing a non-zero 
6 Li/ 7 Li isotopic ratio to a 2 a confidence level. In 
order to account for the microscopic diffusion pro- 
cesses during MS let us recall that our initial 6 Li 
has been increased by 0.2 dex with respect to the 
0.93 dex of the Cayrel et al. (1999) measurement 
in HD84937. Our modelling of pre MS 6 Li history 
suggests features whose consequences are clearly 
not observed in the present halo stars. First the 
depletion increases for decreasing effective tem- 
perature even in the [6500,6000] K range (figure 
3). This trend is not present in the observations, 
which also systematically suggest a higher abun- 
dance than predicted. For instance the metallicity 
of HD106038 (T eff = 5905K) is [Fe/H] =-1.35 dex 
but this star is almost 2 dex above the T e g — 6 Li 
track for [Fe/H]=-1 dex at its effective tempera- 
ture. Also because the outer convection zones ex- 
tend closer to the 6 Li than the 7 Li burning regions 
the depletion appears to depend on the metallicity 
between [Fe/H]=-2 and -1 dex. On the contrary, 

4 In this respect we remark that the evaluation of element 
opacities Table 3 of Piau & Turck-Chieze (2002) is incorrect 
in the very cases of hydrogen and helium. 



observed 6 Li abundances appear independent of 
metallicity : in figure 3 the cases of HD68284 
(coolest object) and HD160617 (object closest to 
6000 K) are striking. Although they have metal- 
licities of -0.59 dex and -1.76 dex respectively they 
display similar 6 Li fractions. The difference we ex- 
pect for a 1 dex variation of [Fe/H] is over 1.5 dex 
if T e ff < 6000K. As a matter of fact the data sug- 
gest a 6 Li plateau with respect to [Fe/H] (Asplund 
2005a; 2005b). We note that there also seems to 
be a plateau in the T e g — 6 Li relation. As the de- 
pletion rates are too large when compared to this 
plateau, the situation of 6 Li in population II stars 
is very similar to the situation of 7 Li in population 
I stars regarding pre MS: the actual computations 
predict an overdepletion (PTC02). Interestingly 
this discrepancy lingers despite very different com- 
positions and the improvements in the physical in- 
puts to the code. 

5. Lithium isotopes during the main se- 
quence 

Pre main sequence computations suggest that 
the 7 Li abundance pattern on the Spite plateau 
or below T e ff = 5500K was not achieved through 
early evolution but during the MS. Furthermore 
6 Li shows no systematic pre MS depiction at least 
down to ~ T c ff = 5900K. In this section we study 
the effect of the MS tachocline mixing (Spiegel & 
Zahn 1992) on both lithium isotopes. Pre MS de- 
pletion has been systematically cancelled. Thus 
we assume that all the stars exhibit on their sur- 
face their formation lithium on the ZAMS. This 
assumption originates from the present computa- 
tions indicating that 6 Li pre MS depletion is over- 
estimated. In order to study correctly the lithium 
isotopes history during MS one should not include 
pre MS effects that are not fully understood at the 
moment. 
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Table 4: Opacity and impact of the main metals contributing to the opacity (0, Ne, Mg, Si, S and Fe) at 
the BCZ. The opacities are computed at the peaks of BCZ temperature (i.e maximum lithium depiction 
rate) in three objects that end up with similar effective temperatures on the ZAMS. The ZAMS effective 
temperatures are 4975 K, 4991 K and 4982 K for the [Fe/H]=-2 dcx, -1 dcx and -0.7 dex respectively. This 
temperature corresponds to a moderate 7 Li depletion during the pre-MS (see fig 2). 





[Fe/H]=-2 


[Fe/H]=-1 


[Fe/H]=-0.7 




M* = 0.58M Q 


M* = 0.62M G 


M* = 0.70M Q 


Pre MS maximal BCZ temperature 


3.44 10 6 K 


3.45 10 6 K 


3.79 10 6 K 


& corresponding density 


3.98g.cm~ 3 


3.95g.cm~ 3 


2.66g.cur 3 


Opacity (k) 


5.37 cm 2 .g _1 


8.95 cm^g- 1 


7.62 crnlg^ 1 


Contribution of 








Nc Mg Si S Fe to k 


7.9 % 


43 % 


61 % 


Contribution of 








C N Na Al P CI Ar 








K Ca Ti Cr Mn Ni to k 


0.4 % 


2.6 % 


3.9 % 



5.1. The 7 Li plateau 

When coupled to BBN computations, two in- 
dependent observational facts allow constraints on 
the primordial 7 Li/ 1 H ratio : the deuterium abun- 
dance measurements in quasar absorption systems 
(Buries 2002 and references therein) and the cos- 
mic microwave background anisotropics probed by 
the Wilkinson Microwave Anisotropy Probe (Cy- 
burt et al. 2003; Coc et al. 2004). These con- 
straints both point towards a similar 7 Li fraction 
~ 3 times higher than what is observed in halo 
dwarfs of effective temperature larger than 5500 
K (see figure 2). The reason for this discrep- 
ancy does not currently seem to originate in uncer- 
tainties on the BBN cross sections (Cyburt et al. 
2004). The lithium depletion during halo dwarfs 
evolution therefore remains a plausible explana- 
tion. Besides this, if one considers population I 
stars' histories, a change in the surface abundances 
is expected for halo dwarfs. Helioseismic mea- 
surements strongly support microscopic diffusion 
in the Sun (Bahcall, Pinsonneault & Wasserburg 
1995) i.e. on a shorter timescale and for a deeper 
convection zone than experienced in Li-plateau ob- 
jects. Moreover, a non standard mechanism is nec- 
essary to explain the time dependent 7 Li depletion 
observed in open-clusters G type stars as well as 
the actual solar photospheric lithium abundance. 



There is no physical reason why such a process 
should not be at work in population II stars. The 
interplay of this process with microscopic diffusion 
sets the global lithium history. 

Several mechanisms have been proposed to ex- 
plain the 7 Li depletion in Galactic disk stars and 
were also applied to halo stars : mass loss (Vau- 
clair & Charbonnel 1995), mixing by internal 
waves (Montalban & Schatzman 2000) and rota- 
tion effects. These latter mechanisms have been 
particularly explored notably as during population 
I MS lifetime both depletion and rotation rates are 
observed to decrease. Recently Theado & Vauclair 
(2001) considered the effects of meridional circula- 
tion coupled to molecular weight gradient effects. 
Their models were able to produce a nearly con- 
stant 7 Li depletion for the stars on the plateau 
without ad hoc adjustment of their parameters. 
Pinsonneault et al. (2002) used a rotationaly in- 
duced mixing model calibrated on the Sun. They 
found that the scatter observed on the plateau and 
the small fraction of stars below the plateau were 
compatible with a mild depletion from 0.1 to 0.2 
dcx. 

Whatever the process responsible for the Spite 
plateau lithium evolution, this process has to com- 
ply with two challenges. First the plateau is re- 
markably flat with T e ff : stars with quite different 
outer convection zones experience similar deple- 
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6600 6400 6200 6000 5800 5600 5400 
Effective temperature (K) 

Fig. 3.— 6 Li pre MS depletion in [Fe/H]=-2 dcx 
(solid line) and [Fe/H]=-1 dex (dashed line). Each 
curve has been shifted upwards by 0.2 dex in order 
to account for microscopic diffusion effects. Dia- 
monds: data in the [Fe/H]= -1 to -2 dex range 
from Asplund et al. (2005a) and M. Asplund pri- 
vate communication; Stars: data from Nissen et 
al. (1999) [Fc/H] w -0.6 dcx. 

tion : Ryan et al. (1996) report a ~ 0.04 dex 7 Li 
variation every hundred Kelvin. Second the scat- 
ter on the plateau is extremly small for field stars. 
Spite et al. (1996) find ar u ~ 0.06 - 0.08 dex, 
Bonifacio & Molaro (1997) find a? u ~ 0.07 dcx, 
while Melendez & Ramirez (2004) find cr7 Li ~ 0.06 
dex and Ryan et al. (1999) find <77 L ; < 0.02 dex. 
Given the observational errors the scatter is in any 
case considered consistent with zero by the authors 
and clearly lies well below 0.1 dex 5 . 

We compute the MS evolution of 7 Li for 
[Fe/H]=-2 dcx models (composition described in 
detail in tables 1 & 2). The rotation law is that 
of case A : K = K Q , and Td = 3 Myr (see §3). 
We stress however that a different early rotation 
history has no impact on the MS tachoclinc mix- 
ing and therefore on the surface lithium evolution. 
The reason stems from the swiftly converging ro- 
tational velocities after a few megayears combined 
with the few megayear the tachocline mixing re- 
quires for its onset. In our approach the sur- 
face abundances arc not sensitive to the ZAMS 
angular momentum content. In spite of a wide 

5 The scatter observed in globular clusters such as M92 is 
higher (Boesgaard et al. 1998) and probably induced by the 
interactions resulting from the denser stellar environment. 



range of possible early rotation histories the late 
MS 7 Li contents will show no resulting scatter. 
Even if we exceptionally consider the effect of the 
tachocline mixing to start as early as 50 Myr, 
differences resulting from early rotation differ- 
ences (i.e. the different initial angular momentum 
amounts) would not be significant. In this situ- 
ation the difference in 7 Li fraction we compute 
between the case of a slow and a fast ZAMS ro- 
tation (case A and B respectively on figure 1) is 
0.04 dex at T e ff = 5270 K and gets smaller above 
this temperature. Figures 4 and 5 present our re- 
sults along with observations. The observational 
data have been taken from Ryan et al. (1996), 
Ryan, Norris & Beers (1999), Ryan et al. (2001), 
Melendez & Ramirez (2004) and Charbonncl & 
Primas (2005). In order to diminish the possi- 
ble 7 Li trends with [Fe/H] we selected stars with 
metallicity between -2.5 and -1.5 for the plateau 
objects (T eff > 5500K). Below T cff = 5500K we 
use the data from Ryan & Deliyannis (1998) and 
Thorburn (1994). To keep a reasonable number 
of stars in this region we considered a less restric- 
tive criterium on the metallicity and extended the 
sample to objects having —3 < [Fe/H] < — 1. We 
carefully excluded the subgiants from the sample 
using the recent work of Charbonnel and Primas 
(2005). This work takes advantage of the HIP- 
PARC OS parallaxes in determining the stellar 
evolutionary status. All the objects considered 
in figures 4 and 5 are dwarfs or turn-off stars 
and comply with our prescriptions on metallic- 
ity. However the four objects having upper limits 
in lithium detection and T c ff above 5900 K have 
-1.66 dex < [Fe/H] < -0.88 dex. Only one of 
them is a confirmed dwarf (see table 5). Never- 
theless because of their high effective temperature 
these stars, if subdwarfs, have not experienced any 
significant additive extension of their outer con- 
vection zone since the turn-off. Thus their surface 
lithium fraction is not expected to have changed 
since the main sequence. The computations have 
been made with both tachocline and microscopic 
diffusion or only microscopic diffusion. In any 
case a 7 Li plateau is predicted above 5500 K. This 
plateau lies between 0.2 and 0.4 dex below the 
initial abundances. In the microscopic diffusion 
models the abundances moreover decrease by 0.2 
dex between 5700 K and the hotter end of the 
plateau near 6400 K. On the plateau the main 
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effect of the tachocline mixing is to correct this ef- 
fective temperature dependence : only a small de- 
crease of less than 0.1 dex in 7 Li remains towards 
the hot edge of the plateau. These computations 
are in agreement with the predictions of Richard 
ct al. (2005): on the one hand the microscopic 
diffusion plays the major role in determining the 
level of the plateau, on the other hand the shape 
of the plateau requires that a non standard mixing 
process occurs in the radiation zone. For the sake 
of clarity the tachocline mixing models built with 
the highest rotation velocity of the radiation zone 
(case D) do not appear on figure 4. Above 6000 K 
these models do not show differences in 7 Li exceed- 
ing 0.1 dex with the predictions of case A rotation 
models. In the framework of the tachocline mix- 
ing the most rapid rotation rates inferred from 
HB fast rotators are thus unable to explain the 
lithium poor stars near the turn off. We shall re- 
turn to this question in the next chapter. Below 
5700 K the tachocline mixing models provide an 
explanation of the observed 7 Li pattern. On the 
contrary the pure microscopic diffusion models do 
not predict any 7 Li depletion around 5500-5400 K. 
The tachocline mixing therefore improves the sit- 
uation with respect to pure microscopic diffusion 
models. However there are still two remaining is- 
sues. First the observed plateau still lies 0.2 dex 
below the predictions. Second the abundances be- 
low 5500 K appear largely scattered at a given 
effective temperature. 

A possible explanation of the discrepancy be- 
tween the observed and predicted plateau may 
stem from a slight underestimation of 7 Li. It has 
been suggested recently by Melendez & Ramirez 
(2004) that such an underestimate could be re- 
lated to an incorrect evaluation of the T e ff scale. 
The 7 Li abundance infered by these authors is 
2.37 dex which is in perfect agreement with our 
predictions as shown on figure 5. These observa- 
tional results however still have to be confirmed. 
What is the origin of the scatter in lithium abun- 
dances at a given T e g? on the cooler side of the 
plateau? As we have seen in §4, metallicity ef- 
fects and pre MS structural evolution are unable 
to explain this scatter. If one then investigates 
the effects of a MS rotational mixing process it is 
tempting to relate the lithium scatter to different 
rotation histories. However, several clues coming 
from population I observations lead us to think 
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Effective temperature (K) 

Fig. 4. — T e ff— 7 Li MS relation in tachocline mod- 
els of buoyancy frequency 10 ^iHz (solid line) and 
2 ^iHz (dashed line). The 7 Li depletion patcrn is 
also provided for pure microscopic diffusion mod- 
els (dot-dashed line). The age of the models is 
13 Gyr, they all have the standard composition 
([Fe/H]=-2 dex). The data are the same as in fig- 
ure 2. 

that this is an incorrect interpretation. For in- 
stance, the Hyades solar-like stars all are slow ro- 
tators with vsini < 10km. s _1 . They presumably 
have experienced very different rotation histories 
as is suggested by the observations of younger clus- 
ters. However they show almost no scatter of 
their lithium abundances at a given temperature 
between 5000 and 6000 K (see for instance fig 6 
of Thorburn et al. 1993). Therefore in the case 
of the Hyades it appears artificial to relate dif- 
ferences of initial angular momentum and differ- 
ences of lithium abundances. One could imagine 
that the Hyades solar analogs have kept radiative 
cores rotating with various speeds from their early 
evolution and that the subsequent loss of angular 
momentum will induce different amounts of mix- 
ing and surface lithium destruction. However the 
presence of fast rotating cores and slow rotating 
envelopes after the Hyades age is not supported by 
hclioseismic or early rotation history constraints. 
As we mentioned in chapter 3, solid body rota- 
tion is probably achieved at the Hyades age. If 
this is correct, the tachocline mixing (or any other 
rotation induced mixing) combined with the mi- 
croscopic diffusion will produce similar effects on 
lithium during the subsequent evolution. We sug- 
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Fig. 5. — T e ff— 7 Li MS relation in tachocline mod- 
els of buoyancy frequency 10 /iHz. The 7 Li deple- 
tion pattern is provided for [Fe/H]=-2 dcx (solid 
line), [Fc/H]=-3 dcx (dot-dashed line), [Fe/H]=-2 
with the solar repartition among metals (dashed 
line). The age is 13 Gyr. Note the significant im- 
pact of the metal repartition on the lithium deple- 
tion. The data have been taken considering the 
stars having —2.5 < [Fc/H] < —1.5 of the sam- 
ple of Melendez & Ramirez (2004) observations 
(star symbols) plus Ryan & Deliyannis (1998) and 
Thorburn (1994) for objects below 5500 K in T cff 
(diamonds). 

gest the lithium scatter originates in differences 
of the metal content instead. The relevant quan- 
tity is not [Fe/H] however: as shown on figure 5 
tachocline models with [Fe/H] =-2 dex or [Fe/H]=- 
3 dex exhibit the same depletion pattern. The 
abundance repartition among metals has a much 
stronger influence on lithium. The models keeping 
[Fe/H]— -2 dex but with solar repartition among 
metals are less lithium depleted for any given effec- 
tive temperature than the models assuming a popu- 
lation II metal distribution. Changing the oxygen 
to iron ratio by 0.3 dex affects the 7 Li content more 
at a given effective temperature than changing the 
global metallicity by 1 dex. Figure 6 displays the 
BCZ temperatures as a function of the effective 
temperatures for the three compositions displayed 
on figure 5. The relation is more affected by a 
change in the metal ratios than by a change in the 
total metal fraction. The solar metal repartition 
models exhibit considerably lower BCZ tempera- 
tures at any given effective temperature. The rea- 




6500 6000 5500 5000 4500 

Effective temperature (K) 

Fig. 6. — T c ff — Tbcz relation in the models of 
[Fe/H] =-2 dex, [Fe/H] =-3 dex and [Fe/H] =-2 with 
the solar repartition among metals. Linestyle con- 
ventions are the same as in figure 5. Note the huge 
impact of metal repartition on the relation. 

son for this behavior is twofold. The iron tunes 
the opacity near the top part of the convection 
zone whereas the oxygen tunes the opacity near 
the base of the convection zone. Figure 7 illus- 
trates these opposite behaviors for a model where 
7 Li depletion becomes significant (T e g = 5270K). 
An increase of the iron fraction tends to increase 
the thermal gradient V = in the superadia- 
batic region and thus tends to increase the entropy 
of the deep convection zone. In turn this decreases 
the depth of the convection zone. If the iron con- 
tent is fixed so is the entropy of the convection 
zone and the position of the BCZ then becomes 
very sensitive to the opacity in this region. The 
oxygen being the main opacity contributor near 
the BCZ, a diminution of the oxygen content at 
fixed [Fe/H] will induce a drop of BCZ tempera- 
ture as illustrated on figure 6. The oxygen to iron 
fraction drops by a factor 2 from the population 
II to the population I metals repartition, this in 
turn translates into a smaller 7 Li depletion (figure 
5). 

Because of the combined opacity roles of oxy- 
gen and iron we predict a correlation between 7 Li 
abundance and [Fe/OJ. We expect this correlation 
only below 5500 K and moreover for objects hav- 
ing [Fc/H] above ~-2 dex. If the T e g condition 
is not fulfilled the BCZ is not close enough to the 
regions where 7 Li nuclear burning can occur. Al- 
ternatively if the [Fe/H] condition is not fulfilled 
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Fig. 7. — Left panel: opacity contributions in % for iron (solid line) and CNO (dashed line) as function 
of the pressure in the upper convection zone. The computations are for the 0.62M© (T e ff = 5260K) at 13 
Gyr with [Fe/H]=-2 dex and halo repartition among metals. The superadiabatic region extends down to 
log P=7. Right panel : opacity contributions for iron and oxygen as a function of the pressure in the lower 
convection zone. The linestyles are similar to the upper panel. The vertical line shows the BCZ. 



the metals do not have enough impact on the stel- 
lar structure to induce differences in the history of 
7 Li. In this respect the anti-correlation between 
[Fe/O] and 7 Li observed in NGC6752 by Pasquini 
et al. (2005) does not contradict our expectations 
because it concerns turn-off stars (T e ff > 5900K) . 
Following these authors we think that the chem- 
ical properties of these objects are related to the 
material that formed them. The case of a recently 
studied star HE1327-2326 (Frebel et al. 2005) also 
drew our attention. With an effective tempera- 
ture of 6180 ±80K, HE1327-2326 should lie on the 
lithium plateau. It is currently the lowest metal- 
licity star known: [Fe/H]=-5.4 dex. However it 
is found to be anomalously lithium poor, the up- 
per limit being [ 7 Li] = 1.6 dex. The metal distribu- 
tion in this object differs significantly from other 
population II dwarfs : e.g. [O/Fe] could be as 
high as 3.7 dex for this star and [C/Fe]=3.9 dex. 
Considering the composition observed in this star 6 
we have tested the tachocline mixing model for 
T ff = 6200 K. Despite a large oxygen fraction, 
we found no significant depletion of 7 Li at 13 Gyr. 
Because of the very small metallicity the effects 
of the metals are quite moderate on HE 1327-2326 



The helium mass fraction was assumed 0.2479 while the 
oxygen and all the metals not mentioned in Frebel et al. 
(2005) were set to [X/Fe]=3.7 dex. 



even though [O/Fe] is huge. This illustrate the 
[Fe/H] condition we mentionned just before. In 
this respect it is interesting to note that the heavy 
elements mass fraction and [O/H] ratio we com- 
pute for HE1327-2326 are Z = 3.1 10~ 4 and -1.8 
respectively which is comparable to our standard 
composition case where Z = 3.0 10~ 4 and [0/H]=- 
1.7 dex (see table 1). 

5.2. Below the 7 Li plateau 

About 7% of the halo turn-off stars are lithium 
poor with abundances below 1.7 dex. Some of 
these objects show peculiar chemical abundance 
patterns while some others cannot be distin- 
guished from the other 'normal' lithium abun- 
dance stars. However the majority of these stars 
clearly rotates faster than their Li-normal counter 
parts (Ryan & Elliott 2005). At first sight this 
point should once again suggest a connection be- 
tween lithium and rotation effects. However we 
do not think that these stars are lithium poor 
because they reached the ZAMS with more an- 
gular momentum. The distribution of rotational 
periods in ZAMS open-clusters such as IC2602 
and IC2391 do not show a small fraction of stars, 
say between 5 and 10 %, having periods clearly 
shorter than the bulk of other comparable objects. 
It rather seems that rotation periods are uniformly 
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distributed on the 0.2-10 days range (Barnes et al. 
1999). If the ZAMS rotation distribution was re- 
sponsible for lithium poor halo dwarfs at turn-off, 
one would expect more stars below the plateau 
and a continuous distribution of them as func- 
tion of the lithium fraction. This is not the case: 
hot lithium poor stars are rare and there is an 
observational gap of objects between [ 7 Li] « 2.1 
dex and the upper limit of 1.7 dex. Furthermore 
the Li-poor objects generally belong to binary 
systems. Indeed three out of four Li-poor stars 
of the Ryan et al. (2002) sample are confirmed 
binaries, the secondary component being presum- 
ably a compact object. This leads the authors 
to assume the higher rotation rates were achieved 
through angular momentum accretion accompa- 
nying a mass transfcrt from the former primary 
component. Several processes could account for 
the low lithium fraction : the variation of rota- 
tion rate and angular momentum of the accreting 
object may have induced enough internal mixing. 
Alternatively the accreted matter may have been 
initially lithium poor and induced a subsequent 
dilution. 

Our purpose here is to briefly check these two 
possibilities in the case where one considers the 
tachoclinc diffusion as the rotation induced mixing 
process. Table 5 shows the Ryan et al. (2002) data 
for their Li-poor stars while the table 6 provides 
the general characteristics of our models computed 
for [Fe/H]=-1 dex, more or less the metallicity of 
the observed objects. Out of the four objects of 
table 5, CD-31 19466 might be peculiar: it is nei- 
ther a confirmed fast rotator nor a spectroscopic 
binary and is further away from turn-off than the 
other stars. On this basis we exclude it from the 
subsequent discussion. Table 6 shows that the 
masses of the observed Li-poor stars probably lie 
between 0.8 and 0.86M Q . Masses above 0.9 M Q 
are excluded unless the stars are evolving towards 
the subgiant stage and already have at least lost 
s=s 160 K since the turn-off. Moreover such high 
masses would mean that the Li-poor objects have 
half the age of the typical halo stars. Masses below 
0.8 M also seem excluded because such objects 
should not reach the high observed effective tem- 
peratures. 

Now that the plausible mass range for lithium 
poor objects is established let us consider the effect 
of the tachocline mixing. Ryan et al. (2002) es- 



timate a typical mass increase of 10 2 M . Thus 
the masses of the lithium poor progenitors were 
probably well above 0.75 M (T off = 5700K on 
the ZAMS). Prior to the mass transfert these stars 
were hot enough to be on the plateau and to have 
a 'normal' 7 Li fraction. The depletion impact of 
the tachoclinc mixing on a 0.8 M© star is negligi- 
ble even for a high angular rotation velocity. We 
built a 0.8 M model with constant rotation rate 
1.27 lO^rad.s -1 on the MS. This rate is chosen 
such that at 13 Gyr the surface equatorial veloc- 
ity is 10 km.s -1 which is the order of magnitude of 
the velocity suggested by the vsini measurements 
(table 5). However this model is extreme in terms 
of rotation and tachocline mixing because it keeps 
a high angular velocity all along its evolution. We 
find that the final 7 Li of this model is 2.32 dex. 
Despite the high rotation rate no depletion occurs 
because the region of nuclear lithium burning lies 
too deep below the base of the convection zone. 

Let us now consider the effect of lithium dilu- 
tion. As can be seen in table 6 the masses of the 
outer convection zones of the considered objects 
are quite small. If the total mass of the 7 Li object 
is close to 0.8 M then a 3 1O~ 2 M mass trans- 
fert of lithium free material would be necessary 
to dilute the 7 Li below the observed upper lim- 
its and down to [ 7 Li] = 1.61 dex. Alternatively if 
the mass is close to 0.86 M then a 1O~ 2 M ac- 
creted mass would be enough to decrease lithium 
to [ 7 Li] = 1.55 dex. Therefore the dilution mech- 
anism of 7 Li can explain the low lithium surface 
abundances observed in some turn-off halo dwarfs. 
The need for an additive mixing process is required 
only if the accreted mass is below 1O _3 M . Then 
the decrease in surface lithium fraction would be 
less than 0.2 dex even in the lighter case we es- 
timate for the convection zone. We note further- 
more that in this case the process responsible for 
the mixing could not possibly be the tachoclinc 
mixing. Abundance measurements of other ele- 
ments than 7 Li would certainly clarify the issue 
of the origin of near turn-off Li deficient stars. 
For instance the detection of higher nitrogen and 
lower carbon fractions than in the 'normal' halo 
dwarfs would support the accretion scenario from 
a red giant star. The conversion of carbon into 
nitrogen is ascribed to deep mixing and is clearly 
observed from moderately (Smith, Briley & Har- 
beck 2005) to extremely (Spite et al. 2005) metal 
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Table 5: Turn off lithium poor objects. Data from Ryan et al. (2001) and Ryan et al. (2002). These stars 
appear in figure 2 and 4. Charbonnel & Primas (2005) showed that Wolf 550 is a dwarf and G202-65 is a 
subgiant. The two other stars are not mentioned by this work. 



Object 


vsini 


T e ff 


[Fe/H] 


7 Li 




(km.s -1 ) 


(K) 


(dex) 


(dex) 


CD-31 19466 


< 2.2 


5986 


-1.66 


< 1.49 


Wolf 550 


5.5 ±0.6 


6269 


-1.56 


< 1.61 


BD +51 1817 


7.6 ±0.3 


6345 


-0.88 


< 1.64 



G202-65 8.3 ±0.4 6390 -1.32 < 1.67 



Table 6: [Fe/H]=-ldex models characteristics at their maximum MS effective temperature. 



Mass 


T ff 


age 


Mass of the convection zone 


(M ) 


(K) 


(Gyr) 


(M ) 


0.80 


6233 


12.15 


6.5 10~ 3 


0.83 


6326 


10.25 


3.6 10~ 3 


0.86 


6426 


8.73 


1.9 10~ 3 


0.90 


6548 


7.10 


710~ 4 



poor giants. The presence of 6 Li in these other- 
wise 7 Li poor objects would also favour the ac- 
cretion scenario because it is difficult to imagine 
how this more fragile isotope could survive if 7 Li 
experienced enhanced proton capture. 

5.3. Predictions on 6 Li 

The region of nuclear destruction of 6 Li is close 
to the lower limit of the outer convection zone 
in halo dwarf stars. In this respect 6 Li pro- 
vides tighter constraints than 7 Li on these ob- 
jects. Figure 8 is the equivalent for 6 Li to fig- 
ure 4 for 7 Li. The data are those mcntionned 
in §4. We stress that not all the stars are con- 
firmed as dwarfs by the recent analysis of Char- 
bonnel & Primas (2005). However the subgiants 
of the sample are all hot enough that no significant 
dilution or nuclear destruction of 6 Li is expected 
from deeper convective zones. Two tachoclinc dif- 
fusion models are displayed for [Fe/H] —-2 and halo 
star metal repartition and for [Fc/H]=-2 and solar 
metal repartition. Above 6000 K the agreement 
with the data is good for the two sets of mod- 
els. Because the turbulence induced by tachoclinc 



mixing is shallow, both models brake diffusion of 
heavy elements around the turn-off without induc- 
ing significant 6 Li destruction here. This is an 
improvement with respect to the computations of 
Richard et al. (2005). Nonetheless below 6000 
K the tachocline models for the composition ex- 
pected for halo stars predicts too strong a deple- 
tion. The presence of 6 Li in stars as cool as 5900 
K appears problematic. However no hasty conclu- 
sion should be drawn because of the complexity 
of the situation. First, non detection stars are not 
represented on the plot. We recall that for the As- 
plund et al. (2005a) sample only 9 out of 24 stars 
yielded to an unambiguous 6 Li detection. Second 
because of the still rare 6 Li detections, the data 
of stars having —2 < [Fc/H] < — 1 dex are com- 
pared to the models. For the sake of simplicity 
we consider a unique 6 Li initial fraction (« 1.1 
dex after correction for diffusion effects) . However 
cosmic ray activity predicts a significant increase 
of 6 Li over the actual metallicity range (Rollindc 
et al. 2005). Finally we expect that the metal 
repartition effects will change the effective tem- 
perature/BCZ temperature relation as aforcmcn- 
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1.5 




6600 6400 6200 6000 5800 5600 5400 
Effective temperature (K) 

Fig. 8. — T e ff — 6 Li relation for tachocline diffusion 
models of buoyancy frequency 10^i Hz. Solid line: 
[Fe/H]=-2 dex and halo repartition among metals. 
Dashed line [Fe/H]=-2 dex and solar repartition 
among metals. The data and error bars are from 
Asplund et al. (2005a) (diamonds) and Nissen et 
al. (1999) (star symbols). 

tioned. In the same manner as in the preceding 
chapter we considered the tachocline mixing im- 
pact when [Fe/H]=-2 dex and the metal reparti- 
tion is solar. As can be seen on figure 8 the agree- 
ment to the data is then quite improved. This 
suggest there should be a [Fe/O] to 6 Li correla- 
tion for halo dwarfs around the turn-off. 

6. Summary 

We examined in detail the evolution of both 
lithium isotopes in halo stars of effective temper- 
atures in the 4700 to 6500 K range. Their de- 
pletion history has been followed during the pre 
main sequence and the main sequence up to the 
age of 13 Gyr and compared with observations. 
The 9 Be history was also followed but no sig- 
nificant evolution of this species is predicted be- 
cause of its higher temperature for proton cap- 
ture compared to lithium. Our models include the 
most recent physics (§2) in terms of equation of 
state and opacities (OPAL tables), nuclear reac- 
tion rates (NACRE compilation) and atmospheric 
boundary conditions (NcxtGen atmosphere mod- 
els) . Special attention was given to the metal frac- 
tions impact through opacity effects: models with 
[Fe/H]=-3, -2, -1 and -0.7 dex have been consid- 
ered. In all cases the helium mass fraction and 7 Li 



number fraction were assumed to be primordial: 
Y=0.2479 and [ 7 Li]=2.62 dex respectively (Coc et 
al. 2004). The initial 6 Li number fraction was as- 
sumed to be 1.13 dex. The microscopic diffusion 
has always been taken into account. Moreover we 
have studied the effects of the tachocline turbulent 
diffusion, a rotationally induced mixing that was 
previously demonstrated to explain lithium and 
beryllium history in the Sun and other population 
I stars. The two free parameters of the tachocline 
diffusion process are calibrated on the Sun (§3). 
The rotation history of a 0.7 M model is used to 
account for the tachocline mixing. This rotation 
history was computed in detail, taking into ac- 
count the structural changes and surface angular 
momentum losses (§3). 

During the pre main sequence all the models 
evolve from a fully convective stage to a radiative 
core / convective envelope structure. Meanwhile 
the temperature at the boundary between these 
two regions reaches a maximum. The peak tem- 
perature increases as the mass of the star decreases 
which results in an enhanced lithium depletion on 
the ZAMS for cooler objects. As is well known, the 
observed 7 Li abundance pattern shows a plateau in 
7 Li vs T e ff from the turn-off down to T e ff « 5500 K 
while below this value the depletion gradually in- 
creases with decreasing effective temperature. We 
showed that this pattern cannot be explained by 
pre main sequence depletion regardless of metal- 
licity below [Fe/FI]=-l dex (§4). The predictions 
of depiction are at least 0.7 dex above the obser- 
vations in the effective temperature range between 
5500 K and 5000 K. Thus an additive non stan- 
dard mixing process in the radiation zone is neces- 
sary to explain the present day observations. It is 
noteworthy that, contrary to the population I case, 
the ZAMS predicted T ff — 7 Li relation does not 
depend on the metal content. 6 Li abundance pro- 
vide additional informations about the pre main 
sequence evolution. The expected depletion of this 
isotope depends on the metallicity and clearly ex- 
ceeds what is observed in some halo dwarfs. This 
issue is similar to what we found in earlier work 
for population I star with respect to 7 Li pre main 
sequence depletion (PTC02). For surface temper- 
atures cooler than 6000 K some stars are observed 
on their late main sequence with a 6 Li content 
at least 1.5 dex above pre main sequence predic- 
tions. This indicates that the initial conditions or 
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the early structural changes of population II stars 
are not understood. 

The pre main sequence results show that a non- 
standard main sequence mixing near the base of 
the convection zone is necessary to explain the 7 Li 
fraction below T cfr = 5500 K. The T cff - 7 Li re- 
lation calculated with microscopic diffusion mod- 
els at 13 Gyr also suggest that a turbulent pro- 
cess prevents too efficient microscopic diffusion 
near the turn-off. The tachoclinc mixing provides 
good results for the 7 Li abundance pattern on the 
cool side of the plateau. Moreover as the phe- 
nomenon creates turbulence below the convection 
zone it also prevents increased depletion where the 
microscopic diffusion timescale becomes shorter. 
Consequently the Spite plateau is well reproduced 
(§5). Its average level lies between 2.4 and 2.3 
dex. This is slightly above the bulk of the ob- 
servations, but in perfect agreement with the re- 
cent results of Melendez & Ramirez (2004). We 
showed that the early main sequence rotation his- 
tory has no impact on 7 Li provided that nearly 
solid body rotation is efficiently enforced in low 
mass stars. There are indeed several observa- 
tional clues that the solid body rotation is rapidly 
achieved in low mass stars. Furthermore using 
constraints from horizontal branch stars rotation 
we showed that even if halo dwarfs keep fast rotat- 
ing cores this should not induce increased lithium 
depletion around turn-off. The tachocline mix- 
ing requires a few gigayears to become effective 
whereas solid body rotating stars achieve similar 
rotation rates after w 0.2 Gyr. Consequently the 
small scatter on the 7 Li plateau is not contradic- 
tory with scattered rotational rates on the ZAMS 
even if one considers that the turbulence below the 
convection zone is connected to rotation. In the 
framework of the tachocline mixing the initial an- 
gular momentum and the surface abundances are 
not correlated. The tachocline mixing is robust 
with respect to the new observations. It predicts 
6 Li abundances in agreement with the observa- 
tions above 6000 K. The rare objects exhibiting 
a 7 Li abundance clearly below the plateau and an 
effective temperature above 6000 K cannot be ex- 
plained in the framework of the tachocline process 
despite their abnormaly high rotation rates. We 
show however that because of their light convec- 
tion zones, accretion of lithium free matter below a 
few percent of solar mass could explain the actual 



observations. The detection of 6 Li or enhanced 
nitrogen content in these objects would be strong 
support in favor of such an accretion scenario from 
an evolved star. We suggest that the scatter in 
7 Li for stars cooler than the red edge of the Spite 
plateau is related to differences in the metal repar- 
tition in the stars instead of different rotation his- 
tories. The variation of the metals ratio from pop- 
ulation II repartition to population I repartition 
at fixed T e ff significantly diminishes lithium de- 
pletion. We predict a correlation between [Fe/O] 
and lithium abundance in both isotopes at any 
given effective temperature below 5500 K for 7 Li 
and 6000 K for 6 Li. These correlations should be 
observed in dwarfs stars where the metal reparti- 
tion goes from halo repartition to disk repartition 
(i.e. around [Fe/H]=-1.5 dex). To our knowledge 
such correlations have not yet been investigated. 
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